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ABSTRACT 

During the first half of the universe's age, a heyday of star-formation must have 
occurred because many massive galaxies are in place after that epoch in cosmic history. 
Our observations with the revolutionary Herschel Space Observator}Q reveal vigorous 
optically obscured star-formation in the ultra-massive hosts of many powerful high- 
redshift 3C quasars and radio galaxies. This symbiotic occurrence of star-formation 
and black hole driven activity is in marked contrast to recent results dealing with 
Herschel observations of X-ray selected active galaxies. Three archetypal radio galaxies, 
at redshifts 1.132, 1.575, and 2.474 are presented here, with inferred star-formation rates 
of hundreds of solar masses per year. A series of spectacular coeval AGN/starburst 
events may have formed these ultra-massive galaxies and their massive central black 
holes during their relatively short lifetimes. 

Subject headings: Galaxies: formation — Galaxies: starburst — Galaxies: high-redshift 
— Infrared: galaxies 



1 Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator 
consortia and with important participation from NASA 
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Introduction 



The most massive galaxies known have stellar masses ~ 5 x 10 11 M , and the fact that they 
already exist at ea rly epochs (jFontana et al.l 120061 ) implies that they must have formed rapidly 
(IDaddi et al.l 120051). durin g just a few Giga-years. They are expe cted to host supermas sive black 
holes (IHaring &: Rixll2004l ) which have periods of active accretion (IDe Breuck et al.ll2010l ). A sym- 
biotic occurrence of modest star-formation and accretion acti vity has been inf erred in the high mass 
host galaxies of intermediate redshift X-ray selected AGN (|Lutz et al.ll2010l ). Additional interest 
for details of the sy mbiosis in radio-lo ud AGN is moreover generated by the proposed negative 
AGN feedback (e.g., iBower et al.ll2006l ). Here we quantitatively address star- formation and black 
hole build-up in the most massive galaxies hosting radio-loud AGN, during their time of formation. 

Powerful high-redshift radio galaxies are ideal tracers of the extreme form of starburst-AGN 
symbiosis for several reasons. Firstly, their central black holes are being fed at a high rate, hence are 
rapidly growing. Secondly, their large extended radio sources allow an unambiguous quantification 
of the accretion power, from the huge radio luminosities; moreover, their radio morphological 
properties permit an estimate of the duration of the AGN and possibly coeval star-formation 
episodes. Thirdly, the presen ce of large reservoirs of molecu lar gas, the necessary fuel for star- 
formation, has been inferred ([Solomon fc Vanden Boutll2005l ) from (sub)-millimetre spectroscopy 
in some radio galaxies (but not in all). Quantification of the ongoing star-formation is provided 
through far- infrared (FIR) observations, since dust absorbs the radiation of the young stars and 



re-emits it at FIR w avelengths. From rest-frame submillimetre photometry (Benfor d et al 



199 



1 

Reuland et al.l 12004 ) and mid-infrared spectral energy distributions, SEDs (jSevmour et alj 12008) 
it was suspected that some powerful high redshift radio galaxies might be experiencing phases of 
copious (dust obscured) star-formation. The implications from the SED measurements, however, 
were uncertain: the amount of cool dust reradiating the emission of the newly formed stars was 
not well constrained, as the rest-frame far-infrared regions and the Rayleigh-Jeans tails lack data. 



The Herschel Space Observatory (jPilbratt et al.ll2010l ) provides full measurement of IR-submm 
SEDs at unprecedented sensitivity. This Letter presents breakthrough results on a sample of 
seventy radio-loud, high-redshift AGN, namely quasars and radio galaxies from the 3C and 4C 
catalogues. With Pi.4GHz ~ 10 27 ' 5 W Hz -1 , these AGN contain the most powerful accreting black 
holes. Inspection of our sample measurements shows FIR detection fractions of roughly two in three 
objects, at both 160/xm and 250/im (from the Herschel PACS and SPIRE instruments, respectively). 
Stressing that the non-detected objects are probably just below the sensitivity limit, we focus here 
on three archetypal objects, reporting on their spectacular on-going star-formation. 



3C 368 ( 


Best et al. 


1997, 


1998a 


b; 


Cham 


aers & Chariot 


1990; 


Dioreovski et al. 1987 


), 3C68.2 


(Best et al.lll997 


, 1998a 


Chambers & Chariot 


1990l: Dioreovski et al.lll988h and 3C 257 ( 


Van Breueel et al. 


1998 


), at redshifts z = 1.132, 1.575, and 2.474 respectively, are archetypal protogalaxies with 



well-documented pr operties. 3C 257 is the highest redshift object in the complete 3CR sample 



(jSpinrad et al.lll985l ). The first studies of 3C68.2 and 3C 368 in the mid and late 1980s, and their 
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identification as protogalaxies marked the birth of the K-z-diagram ([Lilly &: Longairlll984f) and th e 
beginning of extensive observations of their massive hosts in a cosmological context (|Spinradl ll986). 
The presence of a luminous QSO - a Type-1 AGN, obscured from direct view behin d a toroidal 
circumnuclear dust configuration - was proven in the case of 3C 257 from X-ray data (jDerrv et al 



20031 ); similarly the two other objects reveal obscured QSOs, in our as yet unpublished Chandra 
data. 



Observations and results 



Photometric observations of 3C 368, 3C 68.2 and 3C 257 were carried out as part of a Herschel 
Guaranteed Time program - see Table 1. The PACS scan-map mode was employed with cross scans 
in the blue (70/im) and the red (160/Lim) bands, yiel ding angu lar resolution of respectively 5" and 
11". Data reduction was performed within HIPE (jOttl l2010l ) following standard procedures for 
deep field observations, including source masking and high-pass filtering. Both scan directions were 
processed individually and then mosaiced to yield the final map. Aperture photometry (including 
appropriate aperture corrections) of the sources at the known radio core positions was carried out 
to measure source flux densities. Photometric uncertainties were determined by measuring the 
flux in several apertures on empty parts of the background. The SPIRE small map mode was 
employed in the 250 /im (18" resolution), 350 /im (25") and 500 /im (36") bands. Data processing 
followed standard procedures for small map data within HIPE. The photometry was performed 
utilizing a SPIRE source extractor also implemented in HIPE ([Savage fc Oliver! 120071 ). The un- 
certainties of the SPIR E photometry are dominated by confusion noise. Combining the Herschel 
data with our existing ([Haas et al.l 120081 ) Spitzer mid-infrared photometry, we were able to trace 
our science targets as well as other objects in the fields over two decades in wavelength. With 
this approach we ensure proper identifications while avoiding source confusion. Supplementing our 
Herschel and Spitzer data with published photometry we obtain the complete optical-IR-submm 
SEDs. To extract physical parameters, we fitted the data with three components representing typ- 
ical constituents of a radio galaxy SED. In practice, we started with the AGN-powered warm dust 
emission by cy cling through a library o f torus models covering a wide range in the relevant pa- 
rameter space ([Honig &; Kishimotdl2010i ). For each torus model (using edge-on inclinations of 60°, 



75°, or 90°) we added a black body in the optical-NIR for emission from the host galaxy starqfl as 
well as a modified black body of variable temperature (grey-body model with (3 fixed at 1.6) in the 
FIR/submm to account for possible excess emission due to star-formation. Figure 1 shows for each 
source the combination of one of the torus models, scaling for all components, and FIR dust tem- 
perature which overall minimizes the chi-square, i.e., the best fitting model. Errors on the derived 
physical parameters were determined from the distribution of these parameter values for the range 



2 3C257 required the addition of an extra black body of ~1300K. Such a component, which is generally identified 
with hot dust close to the sublimation temperature, is often observed in Type-1 AGN but has not been seen in 
Type- 2s so far. Inclusion or exclusion of this hot dust component does not affect the forthcoming conclusions. 
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of model combinations consistent with the data. Since radio galaxy jets have large inclinations, 
radio core emission is relativistically de-boosted r ather than boosted as in quasars. Core radio data 
([Best et al.l 119971 . Il998bl ; IVan Breugel et al.l Il998l ) indeed show that any non-thermal contribution 
to the submm emission is negligible. While the SPIRE beam also encloses the radio source hot 
spots, their contribution is negligible since even the total integrated radio SED underpredicts the 
FIR/submm by a large margin. As Figure 1 shows, substantial cool dust emission is called for by 
the SED grey bodies, with temperatures of 53 K, 36 K and 37 K, for 3C368, 3C68.2 and 3C 257 
respectively. Using standard cosmology (Ho = 70 km s _1 Mpc -1 , = 0.73, fi m = 0.27) we 
compute source intrinsic properties, which are listed in Table 2. The listed LAGN _va hies originate 
from the torus modeling: they represent the (accretion) luminosities required to power the fitted 
torus emission. The LsF-values were determined by integrating the fitted FIR grey bodies between 
8^m and 1000/Um. 



3. Discussion 



In analogy with other studies (IBenford et al.lll998l : lHatziminaoglou et al.ll2010i ) we attribute the 
torus heating, i.e., the warm dust component, to the obscured AGN: their luminosities are ~ 10 46 



(Barthcl 


1989; 


Derrv et al. 


2003; 


Haas et al. 


2004: 


Ogle et al. 



2006). Attributing the cool dust emission to star- formation, we determine luminosities which are 
comp arable to the AG N luminosities (Table 2), all in excess of 10 12 L . Using standard conver- 
sions (jKennicuttl Il998l ). the star-formation rate (SFR) values are 610, 390 and 770 M yr" 1 , for 
3C 368, 3C 68.2 and 3C 257, respectively. Hence, all three high-z radio galaxies are prodigious star- 
formerjfi Four points are noteworthy: (1) when comparing the SFRs with values inferred from 



uv/optical data ([Chambers Charlotlll990l ) we conclude that most of the star- formation in these 
three objects is strongly obscured in that band, (2) as judged from their radio sizes (this point 
will also return later), these are mature objects rather than young AGN in a transition phase from 
dust obscured infrared galaxies to AGN, (3) 3C 68.2 and 3C 257 have cool dust tem peratures in 
the range of those measured for distant sub- millimetre galaxies (jMagnelli et al.ll2012l ). and (4) as 
judged from the detection statistics (~ 70%), many (if not all) of our sample objects must have 
optically obscured star- formation of comparable strength, i.e., SFRs of hundreds of M yr _1 . The 
high star-formation luminosities for these high-z radio galaxies are consistent with results from deep 
field H erschel surveys targeting a mix of radio-quiet AGN types over a larg e redshift range out to 



z ~ 



3 (jHatziminaoglou et al 



z > 4 QSO studies (jLeipski et al 



2010 



Elbaz et al 



2010 



Mullanev et al.ll2012al ) and with results from 



2010j) , but our detection fraction is higher. 



3 For comparison, we performed the same decom position for the host of the powerful low-redshift radio galaxy 
Cygnus A (3C 405) , using published ( Haas et al1[20 04') SED data, finding considerably lower values: Lagn = 8.2 x 10 44 
erg sec" 1 and L SF = 1.5 x 10 11 L (SFR 26 M Q yr" 1 ). 
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Our observations indicate that many distant powerful radio galaxies display star-formation at 
a level comparable to ultra-luminous infrared galaxies (ULIRGs) which is coeval to their black hole 
activity. That conclusion together with the above mentioned results of He rschel studies de aling 
with powerful radio-quiet AGN are in marked contrast with the conclusion ( Page et al. 20121 ) for 
X-ray selected 1 < z < 3 (radio-quiet) AGN in the CDF-N, namely that ultra-luminous starburst 
activity is not seen in the most powerful AGN. While detailed comparison of these conflicting results 
is beyond the scope of this Letter we point out two possible explanations. Firstly, the hosts of the 
CDF-N AGN are of lower stellar mass than the radio source hosts. We stress that distant 3C and 
4C objects (of which only of order 10 2 are known) represent the highest peaks in the galaxy mass 
distribution. In fact, our observations are in broad agreement with the extrapolated cosmologically 
evolving speci fic SFR, sSFR, i n the h osts of luminous AGN (and star-forming galaxies) recently 
established by iMullanev et al.l (|2012al ) . A - puzzling - inverse correlation betwe en AGN X-ray 
strength and host galaxy mass would be required to explain the iPage et al.l (|2012l ) results within 
the just mentioned sSFR behaviour in combination with our findings for radio-loud AGN. Secondly, 
radio-selected AGN - selected using radio lobe luminosity which is an isotropic property - represent 
a cleaner, more complete AGN subsample than do X-ray-s elected AGN: the latter may lack the 
most obscured sources. The X-ray luminosities reported by lPage et al.1 (120121 ) are uncorrected for 
obscuration, who estimate the effect is small. X-ray observations for the 3C sources reported here 
indicate observed luminosities, assuming a standard power-law spectrum with T = 1.9, of about 
3 x 1 43 erg s~\ 



10 43 erg s 1 (Wilkes et al., 2012, in prep.) and 1.5 x 10 44 er g s 



(Pe rry et al 



20031 ) for 3C 68.2, 3C 368 and 3C 257 respectively, comparable with sources in the lPage et al.l (|2012l ) 



sample. It is well known that low signal-to- noise (S/N) X-ray data provide poor estimates of 



obscu ration and so result in luminosities uncertain by 1-2 dex (jCappi et al.l 120061 : I Wilkes et al. 



2005). Indeed, for all three 3C sources reported here (radio galaxies), the X-ray data indicate 
significant obscuration with estimated Nh ~ 10 23-24 , yielding estimated intrinsic, hard-band X-ray 
luminosities of ~ 10 45 erg s" 1 for 3C 68 . 2 and 3C 368 (Wilkes et al., 2012, in prep.) and 9 x 10 44 erg 
s" 1 for 3C 257 (2-10 keV, berrv et all ([2003)). The low S/N of the X-ray data combined with the 
likely presence of additional, unobscured, radio-j et-linked emission (resp onsible for the typically 3x 
brighter X-ray luminosities in radio-loud quasars. IZamorani et al.l ( 198ll )) make it likely that the Nh 
and intrinsic X-ray luminosity for these three 3C sources remain underestimated. Although X-ray 
emission in radio-quiet AGN is less complex and may be somewhat lower luminosity, the stronger 
bias against highly-obscured sources and the difficulty i n detecting abso rption in the fainter X- 
ray sources, which are those most likely to be obscured (|Kim et al.l 120070 . is very likely to result 
in significantly underestimated X-ray luminosities in the Page et al. X-ray-selected sample. The 
fact that a higher incidence of X-ray absorption in those sources with 250//m detections is also 
reported strengthens this possibility and so questions their main conclusion: that 250 /im emission 
is undetected, indicating that SF is quenched, in the highest luminosity, Lx > 10 erg s , active 
galaxies in their sample. Careful examination of the true nature of the high-z X-ray AGN is required 
to back-up the proposed negative feedback; such feedback is definitely not observed in the radio- 
loudest) AGN. While the duration of the vigorous star-formation phase could be longer than the 
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radio-loud activity episode, with the star-formation starting earlier than and possibly eventually 
being quenched by the AGN (analysis of our full sample of young, mature and old radio sources 
will address that issue), based on our observations there is a substantial symbiotic period. 

In contrast to other AGN, radio galaxies - with their extended radio structure at high incli- 
nation - offer the unique possibility of estimating the ag e of the activity episode from the linear 
size of the structure, adopting a typical (|Best et al.lll995l ) source expansion speed of ten to twenty 
per cent of the speed of light. The average duration of such an episode follows from the maximum 
size of radio galaxies at the relevant cosmic epoch. Then 3C68.2, measuring ~190kpc, has a large 
dimension implying its AGN phase has an age of several million years, while the sizes of 3C 368 
and 3C 257 indicate that they are somewhat younger. Focussing on 3C 68.2, and adopting an AGN 
age of 5 Myr, a symbiotic star- formation phase at the inferred SFR of ~500 M yr -1 , assuming it 
to be constant, would yield ~ 2.5 x 10 9 M . The actual fi gure could be larg er if the star-form ation 
commenced before the AGN phase. The host mass of 3C68.2, 5 x 10 11 M (jBest et al.lll998al ). can 
be formed through ~200 such gas accretion phases, for which there is sufficient time given the age 
of the universe at that redshift of 4 Gyr. The ~200 associated AGN phases will last a total time of 
~ 10 9 yr during which a massive central black hole of ~ 1 x 10 9 M must accrete ( Haring fe Rix 
20041 ). The implied average black hole fuel consumption of 1 M yr" 1 yields an energy output 
~ r/.dm/dt.c 2 which for a typical accretion efficiency r] = 0.1 equals 5 x 10 45 erg sec -1 - about a 
factor seven below the observed AGN luminosity of 3C 68.2. An essential element in this discussion 
is the assumption that the growth of the central black hole during AGN episodes is in phase with 
the build-up of its host galaxy. We cannot exclude the possibility of a more energetic star-forming 
phase (for instance at SFR ~ 10 4 M yr -1 ), but such SFRs remain to be observed in high-redshift 
radio galaxies. However, if the star-forming phase lasts about a factor ten longer than the symbi- 
otic phase, then ~20 such events are needed to build-up the massive galaxy. In that case the mass 
consumption required to build its black hole during the ~20 AGN episodes would be in agreement 
with the measured AGN luminosity. Other forms of symbiotic occurrences can also be envisaged, 
including pure AGN phases not linked to coeval starbursts and the direct capture of stars through 
galaxy-galaxy merging. Whereas low redshift ULIRGs are believed to originate from major mergers, 
high redshift ULIR Gs in the form of submm galaxies might have formed thr ough cold gas accretio n 



(jDave et al,ll2010l ): the same has been postulated for high redshift AGN (jBournaud et al.l 12011 



Even at their extreme stellar masse s, our 3C hosts obey the sSFR(z) behaviour for non-AGN and 
X-ray AGN (IMullaney et al.l l2012al ). However, their X-ray luminosities are higher, by orders of 
magnitudes; moreover they are radio-loud. They may well mark the rare, most extr e me exa mples 
of the universal accretion to star- formation rate ratio postulated by IMullaney et al.l (j2012bl ) . Fol- 
lowing these authors, we contend that the full symbiosis whereby a cold gas inflow driven, vigorous 
star-formation episode leads into (and may be quenched by) an energetic AGN phase, offers an 
attractive scenario for the reason of a common fuel supply, building up both black hole and stellar 
mass. 
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4. Conclusions 

Following implications from submm emission, mid-IR spectra and Ly-a emission, we now know 
unambiguously that a substantial fraction of the hosts of high-redshift radio galaxies are dust and 
gas rich, undergoing one of likely many phases of copious star-formation. Energetically speaking 
they are comparable to ultra- luminous infrared galaxies. Our study shows accretion and star- 
formation activity simultaneously at work in the most massive early-epoch galaxies, building these 
galaxies and their central black holes, forcefully and effectfully. It will be important to extend the 
investigations to radio-loud AGN hosted by less massive galaxies. Radio-loud activity phases must 
play an important role in galaxy formation. 

The Herschel spacecraft was designed, built, tested, and launched under a contract to ESA, 
managed by the Herschel/Planck Project team, by an industrial consortium under the overall re- 
sponsibility of the prime contractor Thales Alenia Space (Cannes), and including Astrium (Friedrichs- 
hafen) responsible for the payload module and for system testing at spacecraft level, Thales Alenia 
Space (Turin) responsible for the service module, and Astrium (Toulouse) responsible for the tele- 
scope, with in excess of a hundred subcontractors. HCSS/HSpot/HIPE is a joint development by 
the Herschel Science Ground Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS and SPIRE consortia. We acknowledge the support and interest of 
the full Guaranteed Time Proposal team, of Dutch liaison astronomer Max Avruch, of colleague 
Karina Caputi, and the comments of an expert referee. 
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Table 1. Observations and measured flux densities. Herschel observations using the PACS and SPIRE instruments yielded 
photometric data in five bands. The 1-sigma errors combine the measurement noise and the background noise. 



Object 


RA(J2000) 


DEC(J2000) 


Obs.date 


Obs.band 


Integr.timc 


Flux density (mjy) 


3C 368 


18 h 05 m 06?36 


+ll°01'32'/5 


2011 


March 22 


70/^m 


160scc 


30.0±2.0 








2011 


March 22 


160/im 


160scc 


56.0±3.0 








2011 


March 28 


250/im 


lllsec 


31.1±7.0 








2011 


March 28 


350/im 


lllsec 


20.1±6.6 








2011 


March 28 


500/im 


lllsec 


<21.0 


3C 68.2 


02 h 34 m 23?86 


+31° 34' 17'.' 5 


2011 


July 10 


70/^m 


160scc 


27.0±6.0 








2011 


July 10 


160/im 


160scc 


40.0±7.0 








2011 


July 31 


250/im 


lllsec 


28.7±7.1 








2011 


July 31 


350/im 


lllsec 


29.6±7.4 








2011 


July 31 


500/im 


lllsec 


19.7±7.0 


3C 257 


ll h 23 m 09?17 


+05°30'19"5 


2011 


June 1 


70/tm 


640sec 


11.6±1.0 








2011 


June 1 


160/im 


640scc 


16. Oil. 7 








2010 


Nov 23 


250//m 


296sec 


22.8±4.9 








2010 


Nov 23 


350/im 


296sec 


20.1±5.7 








2010 


Nov 23 


500//m 


296sec 


20.4±6.8 



Xable 2. AGN and star- formation luminosities. The inferred AGN and star-formation luminosities for the three high-redshift radio galaxies are listed, as well 
as their (projected) radio source sizes and host galaxy masses, taken from the literature (with references). The 1-sigma errors were determined from the 

distribution of the ten best model combination sets. 



Object 


Redshift 


Tcool 


Lagn 


Lsf 


Lsf 


SFR 


Radio size 


M host 


Mhost /radio refs. 






(K) 


(erg sec -1 ) 


(erg sec -1 ) 


(10 12 L Q ) 


(M yr" 1 ) 


(kpc) 


(Mq) 




3C 368 


1.132 


53±1 


(7.3±3.5)xl0 45 


(1.4±0.04)xl0 46 


3.5±0.1 


610 


80 


10 11.6 


Best et al. (1998a,b) 


3C 68.2 


1.575 


36±3 


(3.3±1.0)xl0 46 


(8.6±0.8)xl0 45 


2.2±0.2 


390 


190 


10 11.7 


Best et al. ( 1997, 1998a) 


3C 257 


2.474 


37±3 


(4.1±2.0)xl0 46 


(1.7±0.1)xl0 46 


4.4±0.3 


770 


95 


< io 11 - 7 


De Breuck et al. (20101; Van Breueel et al. (1998) 



I 
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Fig. 1. — Infrared-submm spectral energy distributions. The infrared-submm SEDs of 3C368, 
3C68.2, and 3C257 are shown, with the best-fitting, multi-component model superposed: stellar 
radiation (red), torus radiation (blue), and cool dust radiation (yellow). 3C257 requires an extra 
hot dust component (green) - see the main text. The mismatch at 16/um for 3C 68.2 is possibly 
due to luminous PAH-emission. The error bars reflect ±l<r errors. 



